Anomalous reduction of the Lorenz ratio at the quantum critical point in YbAgGe 
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We report measurements of the electrical and thermal transport on the hexagonal heavy-fermion 
metal YbAgGe for temperatures T > 40 mK and in magnetic fields H \\ afc up to 14 T. This 
distorted Kagome-lattice system displays a series of magnetic states and a quantum critical point 
at He = 4.5 T. The Lorenz ratio L{T)/Lq displays a marked reduction only close to He- A T-linear 
contribution below 120 mK, present at all different fields, allows to extrapolate the Lorenz ratio 
towards T = 0. At the critical field this yields L/Lo — 0.92 ± 0.03, suggesting a violation of the 
Wiedemann-Franz law due to strong inelastic scattering. 

PACS numbers: 71.10.HF,71.27.+a 



The realization of quantum critical points (QCPs), i.e. 
continuous transitions at T = driven by the change 
of a non-thermal parameter, has led to the discovery of 
interesting new phenomena like non-Fermi liquid (NFL) 
behavior and unconventional superconductivity in recent 
years [IJH] ■ Heavy fermion (HF) systems consist of un- 
stable local magnetic moments (often from Ce-,Yb- or 
U-atoms) in metallic environment and are prototype ma- 
terials for the investigation of QCPs. Their ground state 
is determined by the competition between the Kondo- 
and Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion [5] . As a result of their sensitivity to small changes in 
pressure, composition or magnetic field, HF compounds 
represent an ideal playground for the controlled realiza- 
tion of novel quantum phases [H [7] . 

Two different scenarios were proposed to describe 
QCPs in HF metals: the spin-density-wave (SDW) and 
the Kondo-breakdown scenario [SHTU]. A clear distinc- 
tion between them concerns the role of Kondo singlets 
when the system is tuned from the paramagnetic side to- 
wards the QCP |TT]. In the former scenario, the Kondo 
singlets survive across the QCP and the magnetic order 
is of SDW type. This restricts the critical fluctuations 
to the vicinity of "hot spots", while the main part of 
the Fermi surface remains intact. In striking contrast, 
the entire Fermi surface is influenced in the latter case 
and the heavy Landau quasiparticles break up due to the 
disintegration of the Kondo entanglement [11] . 

Electrical and thermal transport measurements are 
well established tools to investigate the nature of 
QCPs [T^HHj- The ultra low-temperature thermal con- 
ductivity K of a metallic sample is related to its electri- 
cal conductivity a through the Wiedemann-Franz (WF) 
law K/aT = 7r^fc^/3e^ = Lq. Since the WF law is in- 
dependent of the quasiparticle mass it can be viewed 
as a touchstone to verify whether Landau quasiparti- 
cles survive at the QCP. The WF law is strictly valid 
at T = K, if the following necessary conditions are ful- 
filled: (i) the charge and energy must be carried by the 
same fermions and (ii) scattering is fully elastic 15 . If 
for example in a one-dimensional (ID) material charge 
neutral spinon excitations are present, the Lorenz ra- 



tio L{T)/Lq with L{T) ~ K/aT can strongly exceed 
1 in the zero-temperature limit. In 3D metals, a vio- 
lation of the WF law may be expected at the Kondo- 
breakdown QCP where well-defined heavy quasiparticles 
disintegrate p^ll3| . 

Experimentally, a violation of the WF law was re- 
ported in the two prototype HF metals CeCoIns and 
YbRhaSia [HHIl]- In tetragonal CeCoIng, the WF law 
was found to be violated with L/Lq « 0.8 at the field- 
induced QCP when the current flows perpendicular to the 
layers, while it turns out to hold for the in-plane trans- 
port [151 US] ■ The violation of the WF law results from 
linear extrapolation of the thermal and electrical resistiv- 
ity towards absolute zero, which, however, has been ques- 
tioned subsequently [TTj. For YbRh2Si2, two recent stud- 
ies on the thermal and electrical conductivity near the 
field-tuned QCP have found rather similar experimental 
data obtained for fields either perpendicular [TS] (down 
to 25 mK) or parallel [H] (down to 40 mK) to the tetrag- 
onal c-direction at the respective critical fields. However, 
these studies strongly differ in the extrapolation of the 
thermal resistance at the critical field for T — ^ 0. While 
the former study attributed an observed downturn of the 
thermal resistance below 0.1 K to overdamped magnons, 
the latter one treated the same downturn as signature 
of the quasiparticle formation. Exclusion or inclusion of 
this feature to the T ^ extrapolation leads to either 
a violation {L/Lq « 0.9) [13] or verification [M] of the 
WF law, respectively. It is therefore badly needed to in- 
vestigate the validity of the WF law in further quantum 
critical HF metals. 

In this paper, we have performed electrical resistiv- 
ity and thermal conductivity measurements on the HF 
metal YbAgGe in order to investigate the quantum crit- 
ical behavior and the validity of the WF law. Yb ions 
in YbAgGe with hexagonal ZrNiAl structure form a dis- 
torted Kagome lattice, where the magnetic interactions 
between the first nearest neighbors are fully frustrated. 
In zero magnetic field, the specific heat displays a broad 
maximum at 1 K and subsequent sharp antiferromag- 
netic phase transition at 0.7 K [18]. The highly complex 
phase diagram in magnetic field, as shown in Fig. 1 [19) . 
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FIG. 1: (Color online). Schematic T-H phase diagram of 
YbAgGe for H \\ ab, summarizing earlier measurements [19) . 
Below 4.5 T, commensurate AF order forms in phase a, then 
changes to incommensurate AF order in phase 6, and finally 
returns back to commensurate AF order in c |21l 122] . A sign 
change of the thermal expansion coefficient (a) appears near 
the border between the regions c and d. Electrical resistivity 
p{T) varies as T" with n ~ 1 in region d, 1 < n < 2 in 
e and n = 2 in region / [TS]. The red circle indicates the 
position of the presumed quantum critical end point (QCEP) 
at He = 4.5 T. 



implies the importance of magnetic frustration. Interest- 
ingly, the isostructural CePdAl exhibits an unusual par- 
tial ordering, with 2/3 of Ce-4/ moments ordered and 
the remaining 1/3 in a paramagnetic state, due to frus- 
trated interactions 201 . For YbAgGe, neutron scattering 
studies have revealed that the a and c phases are ordered 
with a commensurate wave vector (1/3, 0, 1/3) while the 
b phase displays an incommensurate ordering vector (0, 
0, 0.325) 21, 22 . The microscopic nature of all the other 
phases is still unknown. 

At the critical field separating the c and d phases, H^ — 
4.5 T, the thermal and the magnetic Griineisen ratios di- 
verge as the temperature is lowered, indicating quantum 
criticality [19l [23] . The positions of a cross-over feature 
in the Hall effect and a sign change of thermal expan- 
sion (a=0) also extrapolate for T — >■ to 4.5 T, at which 
field a metamagnetic-like signature of the magnetization 
has been found [24]. These recent observations led to 
a speculation that a new type of quantum critical end 
point (QCEP) is caused by a spin-flop transition of local 
moments [121 HSj- This is different to the case of itin- 
erant moment metamagnetism as observed in Sr3Ru2 07 
and CeRu2Si2 i25( ;26 . A further field-induced QCP in 
YbAgGe has been proposed at i? ~ 7.2 T, where both the 
specific heat C/T and thermopower S/T display a loga- 
rithmic divergence towards low temperatures |19j . Our 
measurements, discussed below, indicate that the WF 
law is satisfied in this NFL region e of the phase diagram. 
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FIG. 2: (Color online). Temperature dependence of the elec- 
trical resistivity p{T) of YbAgGe along the c-direction at var- 
ious magnetic fields H \\ ab. 



Remarkably, however, we observe a marked suppression 
of Lorenz ratio L{T)/Lq at if = 4.5 T. Extrapolation to 
T — 7> suggests a moderate violation of the WF law near 
this critical field. 

The heat and charge transport of YbAgGe were mea- 
sured parallel to the hexagonal c-axis at various magnetic 
fields applied along the ab plane. Details are provided in 
the Supplemental Material (SM) [27]. Fig. 2 displays 
the temperature dependence of the electrical resistivity 
down to 30 mK at selected magnetic fields (further data 
see SM). The sample used in the present work exhibits 
a residual resistivity po — 43 /j,r2cm at zero magnetic 
field. However, this relatively high value is not solely 
related to structural disorder, since the residual resistiv- 
ity is strongly reduced by magnetic field and amounts 
to po = 19 /iilcm, corresponding to a residual resistiv- 
ity ratio of RRR = 5.9, at 14 T [57]. In zero field, 
the rapid drop of p{T) at 0.65 K confirms the first or- 
der transition which leads to a clear hysteresis in re- 
sistivity and magnetization measurements |241 I28| . In 
low fields, this transition is suppressed to lower temper- 
atures and finally disappears at 77 = 2 T. At 4.5 T, a 
weak shoulder-like anomaly is found around 0.4 K. At 
higher fields, the electrical resistivity can be described 
by p(T) — pq + AT" with the exponent n ~ 1.60 and 
1.58 for H = 7.2 and 9 T, respectively. The NFL re- 
gion thus largely exceeds the proposed QCP at 7.2 T in 
agreement with the phase diagram constructed by earlier 
resistivity measurements [IH] ■ This may be related to the 
negative low-temperature magnetoresistance dp/dH < 
in this field regime, which saturates only beyond 10 T 
(see SM [HI). 

The thermal conductivity k{T)/T of YbAgGe as a 
function of temperature is displayed in Fig. 3, at tem- 
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FIG. 3: (Color online). Temperature dependence of the ther- 
mal conductivity divided by temperature k{T)/T of YbAgGe 
at different magnetic fields. 



peratures T < 0.3 K and different magnetic fields up to 9 
T. Further data are shown in SM [57]. In general, heavy 
quasiparticles, phonons and magnons could transport en- 
ergy in YbAgGe. We plot k{T)/T, since Bosonic contri- 
butions from phonons and magnons die out at T — ^ 0, 
giving rise to a finite intercept due to electrons. Upon 
cooling to low temperatures, k{T)/T steadily declines 
in all fields and no indication of saturation is observed. 
Unlike in YbRhaSis [H [H] and CeCoIng [H], where 
the total thermal conductivity is dominated by the elec- 
tronic transport and thus k{T)/T increases with decreas- 
ing temperature, the Bosonic contributions in YbAgGe 
cannot be neglected. The relatively high residual resis- 
tivity in YbAgGe strongly reduces the electronic thermal 
conductivity. This is also refiected in the convex shape 
of k{T)/T in particular at low fields, which makes it un- 
feasible to estimate the phonon contribution to thermal 
conductivity by power-law fitting with reasonable expo- 
nent. Therefore it is impossible to analyze the tempera- 
ture dependence of the electronic thermal resistivity, as 
done for GeCoIns and YbRh2Si2. By contrast, we focus 
on the total thermal conductivity and the Lorenz ratio at 
very low temperatures. The overall behavior presented in 
Fig. 3, indicates that n{T)/T is gradually suppressed in 
low fields and reaches a minimum at if = 4.5 T. At higher 
fields, for H = 7.2 and 9 T, k{T)/T becomes strongly en- 
hanced, which corresponds to the pronounced reduction 
of the electrical resistivity in this field regime. 

Using the electrical resistivity p{T) and thermal con- 
ductivity k{T) we can calculate L{T) — njaT. The 
temperature dependence of the normalized Lorenz ratio 
L(T)ILq = (k/T) • p/Lo is displayed in Fig. 4. We 
are interested in the zero-temperature limit. For this 
purpose it is sufficient to calculate L(T) without sub- 
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FIG. 4: (Color online). Temperature dependence of the nor- 
malized Lorenz ratio L(T)/Lo = K,p/{LoT) for various differ- 
ent magnetic fields. The solid lines indicate linear fits to the 
data. Error bars at data points represent the statistical scat- 
tering while those at the lines near absolute zero temperature 
indicate the systematic error for the extrapolated Lorenz ratio 
arising from the finite width of the contacts. 



traction of the unknown additional phonon and magnon 
contributions, since both will disappear as T — )■ 0. We 
turn our attention to very low temperatures (T < 0.15 
K) at which the Bosonic contributions gradually die out. 
In our accessible temperature range above 40 mK, the 
Lorenz ratio L(T)/Lo > 1. Again this behavior, which is 
caused by the relatively high po in YbAgGe, is in sharp 
contrast to the observations in the very clean HF sys- 
tems CeCoIns dH HH and YbRhaSia [HI [Hj for which 
L{T)/Lq < 1 at finite temperatures below 1 K. Em- 
pirically, we observe a linear temperature dependence 
of L{T)/Lq for all investigated fields with no trace of 
a downturn to lowest temperatures (see SM for further 
data at additional fields and towards higher tempera- 
tures [H]). Although we cannot justify this linear tem- 
perature dependence, which arises from the inclusion of 
Bosonic contributions to k(T), it nevertheless allows us 
to obtain the K values by simple linear extrapolation. 
At _ff < 3.5 T, the Lorenz ratio is weakly field depen- 
dent, while we observe a drastic reduction of the low- 
temperature value at the critical field He — 4.5 T. Upon 
further increasing the field, it increases again. 

While the depression of L(50 mK) clearly indicates 
anomalous behavior near the QCP, we now turn to the 
extrapolation of the linear L{T) dependence found at all 
different fields (see also SM [27]) towards T — > 0. It re- 
veals a value significantly below 1 only at the critical field 
with L{T -^ 0)/Lo = 0.92 ± 0.03. The observed reduc- 



tion of L(T) at He therefore suggests a violation of the 
WF law. On the other hand, for H = 7.2 and 9 T, there 
is a clear change of the slope of L{T)/Lq and the value 
of L{T — ?► 0)/Lo extrapolates to 1.00 ± 0.03 proving that 
the WF law holds in this extended NFL region of the 
phase diagram. 

The electrical resistivity p{T) of YbAgGe is almost 
temperature independent below 0.1 K in all fields. There- 
fore, the electronic thermal conductivity Ke/T should be 
constant, as observed e.g. in CeNi2Ge2 [30]. According 
to the relation n/aT = L{T), the slope of the Lorcnz 
ratio reflects the magnitude of the Bosonic thermal con- 
ductivity (Kph and Kmag) due to phonons and magnons. 
The remarkable depression of the slope in L{T)/Lo at 
higher fields {H = 7.2 and 9 T) is probably linked to the 
loss of magnons due to the suppression of the AF order 
since phonons are unaffected by the magnetic field. 

At He = 4.5 T, the extrapolated T = value of 
L{T)/Lo is less than 1, suggesting that the thermal trans- 
port ability is significantly depressed by additional scat- 
tering sources. Let us return our attention to the raw 
data of the electrical resistivity p{T) and thermal con- 
ductivity k{T)/T at He in more detail. Compared to 
the data at iJ = 2 T, p(4.5 T) is significantly reduced 
below 0.2 K, leading to an enhancement of <t. However, 
the thermal conductivity k(4.5 T)/T does not show a 
corresponding increase. The disparate field response of 
a and k/T therefore leads to the distinct suppression of 
L{T)/Lo at the critical field 4.5 T compared to all other 
fields. This experimental observation is independent on 
the empirical T — >■ extrapolation and indicates anoma- 
lous behavior near the QCP. 

It is interesting to compare YbAgGe with the itinerant 
metamagnets Sr3Ru207 25j and CeRu2Si2 [31]. For both 
systems thermodynamic indications for quantum critical- 
ity in accordance with the itinerant Hertz-Millis scenario 
have been found [26l |32] , although the approach of the 
QCEP is cut-off by the formation of a nematic phase in 
the former and crossover to Landau Fermi liquid behav- 
ior in the latter case. Most interestingly, both materials 
obey the Wiedemann-Franz law within 5% error at their 
quantum critical regions [53J [33] . The most significant 
difference between the QCEPs in YbAgGe and these two 
materials concerns the nature of the magnetic moments. 
In the above mentioned itinerant systems the quasiparti- 
cles stays intact [2nil3i]- Relatedly the WF law is obeyed. 
For YbAgGe, on the other hand, the metamagnetic be- 
havior is more local in nature, resulting in strong critical 
fluctuations that may lead to a destruction of the quasi- 
particles at the entire Fermi surface. At H ~ 4.5 T, 
previous transport measurements have found an abrupt 
jump in the field-dependence of the Hall coefficient [36] 
and sign change of the isothermal thermopower S{H) at 
low temperatures [37] . These observations indeed suggest 



a significant change of the Fermi surface at this field. 

Another distinct difference of YbAgGe to Sr3Ru207 
and CeRu2Si2 is the geometrical frustration of Yb"^"*" 
magnetic moments on the distorted Kagome lattice. It 
is anticipated that the metamagnetic transition is associ- 
ated with a weakening of the magnetic frustration. Once 
the metamagnetic transition is tuned towards zero tem- 
perature, the long-wavelength critical fluctuations of the 
order parameter start showing quantum mechanical char- 
acter and play an important role. The coexistence of 
magnetic frustration and quantum fluctuations may re- 
inforce each other, giving rise to an increasing proportion 
of inelastic scattering of electrons. Close to the quantum 
critical region, strong inelastic scattering of electrons can 
then lead to inefficient heat transport compared to its 
counterpart charge transport, resulting in the violation 
of the WF law. 

A further field-induced QCP in YbAgGe was proposed 
by Schmiedeshoff et al. near 7.2 T, evidenced by the 
onset of a divergence of the Griineisen ratio [TH]. This 
may account for the finite field region between 7 and 9 T, 
for which NFL behavior has been found. In the T — K 
limit, however, the verification of WF law at H = 7.2 T 
implies that the quasiparticles still survive in this NFL 
region. 

In summary, we have studied the low-temperature elec- 
trical resistivity and thermal conductivity of YbAgGe. 
The temperature dependence of the Lorenz ratio 
L{T)/Lq indicates an anomalous suppression at the crit- 
ical field He = 4.5 T of a metamagnetic quantum criti- 
cal end point. Although we cannot neglect or subtract 
the phonon and magnon contributions to thermal con- 
ductivity, the observed linear temperature dependence 
of L{T)/Lq for all different fields allows to extrapolate 
towards T — >■ 0. Such extrapolation yields a value sig- 
nificantly below 1 {L/Lq = 0.92 ± 0.03) only very close 
to He- This suggests strong inelastic scattering even at 
absolute zero temperature due to the disintegration of 
heavy quasiparticles. We speculate that magnetic frus- 
tration and the more local magnetic character compared 
to itinerant metamagnets like Sr3Ru207 and CeRu2Si2 
is responsible for this behavior. 
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Supplemental Material 



I. EXPERIMENTAL SETUP 

Single crystals of YbAgGe were grown from a AgGe- 
rich ternary solution |S1| . Our sample has a hexagonal 
cross section, with ^ 3 mm length parallel to the c-axis, 
along which the heat and charge transport has been real- 
ized. Contacts were made directly on the sample surfaces 
with silver epoxy and were used for both resistivity and 
thermal conductivity measurements. The contacts are 
metallic with typical resistances of less than 200 mil at 
room temperature, leading to a ratio of the contact to 
the sample resistance of about 3. However, this ratio 
should drop below 1 at low temperatures, because of the 
relatively small residual resistivity ratio of YbAgGe. It 
has been shown for silver epoxy contacts on high-purity 
CeCoIns, that the electron-phonon decoupling leads to a 
downturn of k{T)/T at very low temperatures |S2j . Such 
behavior has been absent in our experiments. For the two 
neighboring voltage contacts, the ratio of each contact 
width over their distance is approximately 6%, which can 
be taken as systematic standard error in our measure- 
ments. A dilution refrigerator with a superconducting 
magnet has been utilized, where the field has been ap- 
plied along the ab plane perpendicular to the electrical 
and heat current, respectively. The electrical resistiv- 
ity and thermal conductivity have been measured using 
home-made probes. For electrical resistivity, the base 
temperature has been 30 mK. For the thermal conduc- 
tivity a standard four-wire steady-state method with two 
Ru02 chip thermometers, calibrated in situ against a ref- 
erence Ru02 thermometer located in a field-compensated 
zone has been utilized. The temperature gradient has 
been set to between 5% and 7% and the lowest accessi- 
ble average temperature in these measurements has been 
40 mK. Magnetic fields were applied along the ab plane 
perpendicular to the electrical and heat current. 



II. ELECTRICAL RESISTIVITY 
MEASUREMENTS 

Electrical resistivity measurements as function of tem- 
perature are shown Fig. 2 of the main part of the paper. 
In Fig. SI, we magnify the temperature dependence of 
the electrical resistivity for all the different investigated 
magnetic fields up to 6 T. 

Fig. S2 shows the isothermal field dependence of resis- 
tivity p{H) for magnetic fields up to 5 T at various low 
temperatures. Two clear features can be seen. The first 
one appears at i? ~ 1 T, the second in the vicinity of 
critical field (~ 4.5 T), in agreement with the previous 
resistivity measurements measured down to 0.4 K |S3j . 
The resistivity anomaly at iJ ~ 4.5 T may be related 
to a metamagnetic crossover between two different mag- 
netically ordered phases |S3| . although it remains con- 
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FIG. SI: Temperature dependence of the electrical resistivity 
p(r) of YbAgGe along the c-direction at various magnetic 
fields H II ab. 
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FIG. S3: Temperature dependence of the thermal conductiv- 
ity as k/T vs. T of YbAgGe along the c-direction at various 
magnetic fields H \\ ab. 



shows a relatively weak field dependence. At high tem- 
peratures, T = 400 mK, p{H) increases slowly with in- 
creasing magnetic field and displays a broad maximum 
around the critical field, in agreement with earlier results 
[S3j . Upon decreasing the temperature, this maximum is 
suppressed and a negative magnetoresistance already at 
4 T develops. In the low-field regime below 1 T, anoma- 
lous behavior is discovered at low temperatures. As can 
be seen in the left inset of Fig. S4, several sudden jumps 
are found in p{H). While we do not know the origin 
of these jumps, the observed hysteresis points at some 
variation of micromagnetic domains. The latter may be 
related to the ordering in the a phase, resulting from 
different equivalent directions of the propagation vector 
k = (1/3, 0, 1/3) |S6| or alternatively to the first-order 
nature of the a-b transition at low temperatures [S7j . 



III. THERMAL CONDUCTIVITY 

MEASUREMENTS 



FIG. S2: (Color online). Field dependence oi p{H) isotherms 
for YbAgGe at constant temperatures between 70 and 400 
mK. The left inset enlarges the low-field region for 200 mK. 
The arrows indicate the direction of the field sweeps. The 
right inset displays the 200 mK data at fields up to 14 T. The 
arrow indicates the critical field He — 4.5 T. 



In the main part of the paper, the thermal conductivity 
results are displayed for a selection of applied magnetic 
fields and at temperatures below 0.3 K. In Fig. S3, the 
full data sets up to 0.95 K are shown, at all different 
applied fields. 



tinuous and broad down to 70 mK. Similar behavior has 
also been observed near the metamagnetic transition in 
CeRu2Si2 [S4], while for Sr3Ru207 the formation of an 
electronic nematic phase leads to much sharper signa- 
tures in the isothermal magnetoresistance |S5] . At in- 
termediate fields (1 T < H < 4 T), the resistivity p{H) 



IV. ANALYSIS OF THE LORENZ RATIO 

Fig. S4 displays the temperature dependence of the 
normalized Lorenz ratio L{T)/Lq — (k/T) ■ p/Lq with 
Lq — (""^/S) • (fcs/e)^. The plot displays all measured 
data at temperatures below 0.25, where data for fields 
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up to 4.5 T are shown in part a and data from 4.5 T 
up to 9 T are shown in part b. The hnes display hnear 
fits from the lowest measured temperature up to a vari- 
able maximal temperature for which T-linear behavior is 
found. A significant deviation of the WF law is observed 
at the critical field of 4.5 T, for which the extrapolation 
reveals L/L^ = 0.92 ± 0.03. In comparison to all other 
investigated fields, the data at the critical field indicate 
a pronounced depression of the Lorenz ratio L{T) near 
50 niK. 



FIG. S4: Temperature dependence of the normalized Lorenz 
ratio L{T)/Lo — hp/{LqT) for all different applied magnetic 
fields. The solid lines indicate linear fits to the data. Error 
bars at data points represent the statistical scattering while 
those at the lines near absolute zero temperature indicate the 
systematic error for the extrapolated Lorenz ratio arising from 
the finite width of the contacts. 



[SI] E. Morosan et al, J. Magn. Magn. Mater. 277, 298 

(2004). 
[S2] M.A. Tanatar, J. Paglione, C. Petrovic and L. Taillefer, 

Science 316, 1320 (2007). 
[S3] S. L. Bud'ko, E. Morosan, P. C. Canfield, Phys. Rev. B 

69, 014415 (2004). 
[S4] R. Daou, C. Bergemann, S. R. Julian, Phys. Rev. Lett. 



96 026401 (2006). 
[S5] S. A. Grigera et al, Science 294, 329 (2001). 
[S6] B. Fak et al, J. Phys.; Gondens. Matter 17, 301 (2005). 
[S7] P.G. Niklowitz, G. Knebel, J. Flouquet, S.L. Bud'ko, 

P.G. Ganfield, Phys. Rev. B 73, 125101 (2006). 



